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00 ■ Abstract 

o 

O I The particular properties of synchronization are discussed for the 

system of coupled van der Pol - Duflfing oscillators. The arrangement 
C , of synchronization tongues and the particular properties of their in- 

^ I ternal structure in the parameter space are revealed. The features of 

' attractors in the phase space and in the Poincare section are consi- 

dered. 

q' 

. 1. Introduction 

d 

The system of two coupled van der Pol oscillators is one of the basic mo- 
dels of nonlinear dynamics demonstrating the phenomenon of mutual syn- 
chronization. There are many papers on this theme because this system 
demonstrates a lot of interesting oscillation regimes and types of behavior, 
C\l I such as synchronous and quasiperiodic regimes, synchronization v^^ith dif- 
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ferent phase conditions [Pikovsky et al., 2001], the "oscillation death" effect 
[Aronson et al., 1990; Pikovsky et al., 2001], global bifurcations [Chakraborty 
\Q . & Rand, 1988], etc. Moreover, the results of this system investigation may be 

O I used in the analysis of different electronic, biological and chemical systems 

[Cohen & Neu, 1979; Minorsky, 1962; Neu, 1979; Pavlidis, 1973; Poliashenko 
et al., 1991]. Interest in this problem does not decrease because of its signi- 
ficance and dynamics variety. In this respect we may point out, for example, 
the paper [Ivanchenko et al., 2004], where the resumptive study of this sys- 
5^ i tem is undertaken within the framework of quasiharmonic approximation. In 

particular, they took into account the possible non-identity of oscillators and 
combined nature of coupling. 

The system of differential equations describing the interaction between 
van der Pol oscillators is of the form 

d X , Q - d/Oc , dx dy « 

-d?-^^-y^Tt^^'^^^y^^^Tt-^^='- 
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Here A is parameter characterizing the excess above the threshold of the 
Andronov - Hopf bifurcation in autonomous oscillators, 6 is the frequency 
mismatch between the second and first oscillators, and /i is the coefficient of 
dissipative coupling. 

Investigation of the system of interacting autooscillatory systems may be 
carried out on different "levels of sophistication": by means of Adler equa- 
tion for oscillators phase difference dynamics [Pikovsky et al., 2001; Rand & 
Holmes, 1980], by means of abridged equations, which are correct in quasi- 
harmonic approximation, when oscillators exceed slightly the threshold of the 
Andronov - Hopf bifurcation [Aronson et al., 1990; Ivanchenko et al., 2004; 
Pikovsky et al., 2001], and, at last, by means of initial Eqs. (1) [Chakraborty 
& Rand ,1988; Pastor-Diaz & Lopez-Fraguas, 1995; Storti & Rand, 1982]. 
The last is the most difficult task due to both more complicated behavior of 
the initial system (1) and the solution dependence on the additional para- 
meter 

When studying the dynamics of excited and coupled systems it is impor- 
tant to take into account anisochronism of small oscillations. If we talk about 
the phase dynamics of individual oscillator, anisochronism is the dependence 
of the rate of phase change on oscillation amplitude. This effect may be con- 
sidered, if we insert additional cubic nonlinearity of Duffing oscillator type 
into Eq. (1). Then the equations for coupled oscillators are 



(2) 



dt dt dt dt' 

The parameter P is responsible for anisochronism of small oscillation^. To 
demonstrate this in an explicit form we give the corresponding system of 
abridged equations for the amplitudes of oscillators R and r and the difference 
of their phases [Pikovsky et al., 2001], which may be derived from (2) in 



^For abridged equations this parameter is not essential and can be eliminated due to 
renormalization [Pikovsky et al., 2001]. 

^In system of Eqs. (1) anisochronism is possible also due to not small values of pa- 
rameter A. But we will call for convenience parameter (3 also the parameter of phase 
nonlinearity or anisochronism. 
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the quasiharmonic approximation: 



dR 

'dt 
dr 

di 

dip 

'dt 



r{\ — ji) — + jiR cos (/?, 

5 + ?,(5{r''-R^)-n{^ + -)sin^. 



R{\ -n)-R^ 



+ jir cos </?, 



(3) 



One can see that the rate of change of oscillators' phase difference tp depends 
on the difference of squared amplitudes R and r even if the coupling is absent. 

There are two reasons why taking anisochronism into account is impor- 
tant. First, with taking this factor into account the system equations may 
be brought to the normal form of Andronov Hopf bifurcation in individual 
oscillator [Pikovsky et al., 2001]. In these regards the van der Pol system 
is partly "degenerated" and only the transition to the van der Pol - Duf- 
fing system removes this degeneracy. At the same time choice of the system 
which can be reduced to the normal form of Andronov - Hopf bifurcation is 
important with regard to generalizations and to application of the results to 
another systems. 

Second, it is known that for the case of synchronization by external force 
the anisochronism leads to new dynamics features. It is more obvious for the 
systems with pulse excitation. In that case the isochronous system leads to 
the map for the phase, which was studied by Ding [1986-1988], Glass et. al. 
[1983], Glass & Sun [1994], Keener & Glass [1984], Ullmann & Caldas [1996], 
Viana & Batista [1998]: 



Here C and T are amplitude and period of influence correspondingly. In the 
case of strong anisochronism it may be reduced to the standard circle map 
[Pikovsky et al., 2001]: 



where O = T{l-\-2>(3\/2) is normalized period accounting for additional phase 
lag. 

These two types of maps give two pictures of synchronization tongues 
different in many respects especially in internal structure of these tongues 
[Ding, 1986-1988; Glass et. al., 1983; Glass k Sun, 1994; Keener & Glass, 
1984; Pikovsky et al., 2001; Ullmann & Caldas, 1996; Viana & Batista, 
1998]. This signalizes that the arrangement of synchronization tongues and 





^„+i = ^„ + Q-3/3Csin^, 







3 



their internal structure will demonstrate essential peculiarities for coupled 
anisochronous systems too. At the same time the system of coupled van der 
Pol - Duffing oscillators (2) practically was not investigated in this context. 
One should emphasize that quasiharmonic approximation [Aronson et al., 
1990; Ivanchenko et al., 2004; Pikovsky et al., 2001] is not efficient in that 
case, since it does not describe higher order synchronization tongues and the 
possibility of period doublings and chaos inside them. In the present paper 
we shall give some results of investigation of system (2). 

2. Structure of the parameters space and illustrations of dynamics 

Sufficiently simple and informative method of studying the parameters 
space structure of nonlinear systems is the method of dynamic regime chart 
construction [Kuznetsov et al., 1997; Kuznetsov et al., 2006]. Within the 
framework of such a method we shall mark the oscillation period of the system 
of coupled oscillators by means of different colors on the parameter plane 
(frequency mismatch 5 - coupling value /x). White color corresponds to the 
chaotic or quasiperiodic motions. Cycles periods were calculated by means of 
the Poincare section method: this is the number of points of intersection of 
the phase trajectory on the attractor and the surface chosen as the Poincare 
section. Only those crossings were taken into account that correspond to the 
trajectories coming to the surface from the one side. 

The system under investigation is characterized by four-dimensional phase 
space {x,x,y,y). Therefore three-dimensional hypersurface that is preset by 
means of some additional condition, e.g., zero velocity of the second oscillator 
y = may serve as the Poincare section. In that case the number n of points 
of intersection of trajectory and the hypersurface was determined. Colors on 
the charts are chosen in accordance with the period n. 

The chart of dynamic regimes obtained in such a way for the van der Pol 
oscillators (1) is given in Fig. 1(a) on the (5, /i)-plane for A = 2.5. The value 
of the control parameter is chosen large, that quasiharmonic approximation 
is not efficient. 

We shall characterize mutual oscillations of oscillators also by means of 
the rotation number w. By analogy with [Kuznetsov et al., 2006; Postnov 
et al., 2005] we find numerically the average return time Ty for the chosen 
Poincare section y = and the average return time Tx for the section of the 
first oscillator defined hj x — 0. Then the rotation number is defined as 
w — Ty/Tx- On the chart in Fig. 1(a) there are shown: main synchronization 
tongue with the frequency ratio of 1/1; area of quasiperiodic regimes with the 
embedded system of higher order synchronization tongues among which the 
tongues with the rotation numbers of 2/3 and 3/5 are the most characteristic; 
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the area of the "oscillation death" effect [Aronson et al., 1990; Pikovsky et 
al., 2001] which corresponds to stabilization of the equilibrium state point at 
the origin due to sufficiently strong dissipative coupling. 




b) 4.00 



oscillation 
death area 






Figure 1: Dynamic regime charts for the system (2) for A = 2.5 in case of 
(a) isochronous (/5 = 0), and (b, c) anisochronous (/? = 1) oscillators. 

Now we pass on to the van der Pol - Duffing model (2). Corresponding 
chart of dynamic regimes and its fragment are shown in Fig. l(b, c) for /3 = 1 
that meet the case of sufficiently great anisochronism. 
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The graph of the rotation number w versus the frequency mismatch 6 is 
shown for this case in Fig. 2(a) for different values of the coupling parameter. 
On this graph one can see characteristic wide "steps" that correspond to 
different rational rotation numbers. When coupling parameter is not great 
{fj, = 0.5) steps are feebly marked. This means that quasiperiodic regimes 
predominate. With growth of the coupling parameter (/i = 1) a number of 
sufficiently wide steps emerge. They correspond (in decreasing order) to the 
rotation number values w = 1/2, 3/5, 2/3, etc. If we continue to increase 
the value of the coupling parameter (/x = 2.3), synchronization regime with 
the rotation number of 1/3 begins to predominate. 




Figure 2: Dependences of the rotation number w upon the frequency mis- 
match 5 for A = 2.5, (3 = 1 and different values of the coupling parameter 
1^. 
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One can see in Fig. 1(b) that the presence of anisochronism leads to the 
displacement of synchronization tongues towards the greater values of fre- 
quency mismatch. Synchronization tongues become so wide that one can see 
the situation of their overlapping, which is characteristic for standard circle 
map. Internal structure of this tongues changes too: period doublings and 
transition to chaotic dynamics inside them are observed. It may be seen 
clearly in enlarged fragment of the dynamic regime chart, shown in Fig. 1(c). 

With the increase of the parameter of anisochronism (3 an onset of new "is- 
lands" of doubled periods takes place inside the synchronization tongues lo- 
cating at greater values of frequency mismatch 5 (see Fig. 3). Synchronization 





Figure 3: Dynamic regime charts for the system (2) for A = 2.5 and (a) 
(3 = 1.7; (b) p = 3. 

tongues change their shape with the further increase of parameter f3. They 
fall into two similar in structure pieces, inside which structures called "cross- 
road area" exist [Carcasses et al., 1991; Kuznetsov et al., 1997]. These struc- 
tures are characteristic for systems with period doublings. Practically, two 
systems of synchronization tongues may be observed in the parameters plane: 
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tops of the first one are disposed along the hne of zero couphng, tops of the 
second one - along the boundary of main synchronization area 1/1, which 
is the line of Neimark - Sacker bifurcation. Chaotic area appears between 
these two systems of synchronization tongues. 

When we analyze chaotic areas it is useful to supplement dynamic regime 
charts with the chart of the first Lyapunov exponent (see Fig. 4). Blue color 
shadings correspond on this chart to the synchronization areas with negative 
Lyapunov exponent. Deeper color corresponds to greater absolute value of 
the exponent. Yellow and orange colors hues correspond to exponent val- 
ues close to zero and designate the areas of quasiperiodic regimes. Due to 
finite accuracy of numerical simulations it is not possible to detect perfectly 
zero value of the first Lyapunov exponent. Therefore areas of quasiperiodic 
regimes have slightly varying shadings. Grey color marks the oscillation death 
area. Red color corresponds to chaotic regimes with positive Lyapunov expo- 
nent. The existence of red color on the chart of the first Lyapunov exponent 
verifies presence of chaotic behavior in the system under investigation. 




am s 30.00 



Figure 4: Chart of the first Lyapunov coefficient for the system of coupled 
van der Pol - Duffing oscillators for A = 2.5, /9 = 1. 

Phase portraits of attractors computed at several characteristic points 
within domains of quasiperiodic regimes are given in Fig. 5 on (x, x) and 
(y, ij) planes for the first and the second oscillators correspondingly. To 
put it more precisely these portraits are projections of the attractor in phase 
space (x, X, y, y) onto corresponding planes (x, x) and {y, y). For small values 
of coupling parameter it may be seen that portraits of attractors look like 
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Figure 5: Phase plane portraits at characteristic points of area of quasiperi- 
odic regimes on the dynamic regime chart for A = 2.5, (3 = 1. 



slightly disturbed limit cycles of individual oscillators. At the same time the 
trajectory never approaches the neighborhood of the point of origin. With the 
increase of the coupling parameter the trajectories of individual oscillators 
become more perturbed and there exist a kind of a threshold value when 
the trajectory of the first oscillator may achieve the neighborhood of the 
origin. Corresponding phase portrait looks like the area completely "filled" 
with trajectories. Evidently, this situation corresponds to the case of ill- 
determined phase of this oscillator. When approaching the oscillation death 
area an essential downsizing of attractors occurs (see scales in coordinate 
axis). The orbit ceases again to achieve neighborhood of the origin, but 
structure of phase portraits and shapes of attractors are different from the 
case of small coupling value. 

Three-dimensional illustrations of the dynamics in the "Poincare sections" 
are given in Fig. 6. They illustrate dynamics of the system of coupled oscil- 
lators at the points in the parameters space inside the domains of quasiperi- 
odic regimes and near the threshold of chaos. We should remind that the 
Poincare section was chosen as the hypersurface y = in four-dimensional 
phase space {x, x, y, y). Hence the dynamics of the Poincare return map may 
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be represented by the orbit in three-dimensional space (x, x, y) as is shown 
in Fig. 6. One can see that phase portraits in such Poincare section are in- 
variant curves, which illustrates quasiperiodic character of dynamics. When 
approaching the chaotic area the curves undergo deformations and break 
down via loss of smoothness. 




Figure 6: Three-dimensional "Poincare sections" for the system (2) for 
A = 2.5, (3 = 1. 

Fig. 7 illustrates the dependence of synchronization picture on the param- 
eter of anisochronism (3. There are presented dynamic regime charts on the 
parameter plane (frequency mismatch 6 - parameter of anisochronism /3) at 
fixed coupling parameter values. One can see the "fan pattern" of high or- 
der synchronization tongues spreading away from the main synchronization 
tongue. With the increase of the coupling parameter these tongues widen 
and onset of "islands" of doubled periods takes place inside them. With 
further increase of coupling parameter "crossroad area" structures emerge 
inside synchronization tongues. Phase portraits computed at characteristic 
points of the dynamic regime chart are given in Fig. 7(c). These portraits 
illustrate cascade of period doubling and appearance of chaos. 
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Figure 7: Dynamic regime charts for the system (2) for A = 2.5 and (a) 
H = 1] {h) ji = 1.5; (c) = 2, and phase plane portraits computed at typical 
areas. 

3. Conclusions 

The analysis of synchronization in the system of coupled anisochronous 
van der Pol - Duffing oscillators reveals an essential role that the anisochro- 
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nism plays in the formation of the picture of emerging regimes. Namely, 
in the present paper we have shown that synchronization tongues become 
essentially wider and the situation of their overlapping occurs. Islands of 
doubled periods, "crossroad area" structures and chaotic areas appear inside 
the synchronization tongues. In three-dimensional Poincare section we have 
observed invariant curves which may break down via the loss of smoothness 
with the formation of chaotic attractors. 
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